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Abstract 

We describe computer experiments suggesting that there is an infinite 
family Cl of Riemann zeta cycles A of each size L m 1, 2, 3, ... such that: 
(1) for any non-trivial Riemann zero p and any A G Cl, each A G A is 
the center of a spiral s Pj \ such that (a) the spiral is nearly logarithmic in 
a sense that we will make precise, (b) the point set IJaga s pA con t a i ns a 
subset B P: a of the backward orbit of p containing in turn a pre-image of p 
under each zeta iterate, that (except for p itself) is preserved by zeta, and 
(c) the members of B Pj a are distributed nearly uniformly on s Pj \, in a sense 
that we will make precise. (2) The forward zeta-orbit of any element of 
the backward orbit of p eventually lands on such a B P: a. If the arc length 
(7 Pj \ between A and p of a spiral s Pj \ were finite, and we had determined 
A, cr p , a and (modulo translation) s P: a, then we would have determined the 
position of p. We collect evidence that the s Pj \ are logarithmic in the 
limit as z G s p , a — >■ A, so that the condition cr Pj \ < oo would be plausible. 
We make conjectures on the statistics of the s p> a, and on relationships 
among the A, the backward orbits, and the basin of attraction under zeta- 
iteration of the attracting zeta fixed point 0 ~ —.295905. We do not 
attempt to determine the A, or s P: \ as a function of p and A. We do not 
prove any theorems. Our experiments were done with Mathematica and 
spot-checked with Sage. 


1 INTRODUCTION 
1.1 Overview. 

In this article, we use computer experiments to study the backward Riemann 
zeta orbits of the non-trivial Riemann zeros p, which appear to lie on spirals s p ,\ 
centered upon repelling fixed points A of iterates of the Riemann zeta function. 
We would like to determine these spirals completely (modulo translation). We 
do not get that far in this article. The reason we would like to do so is just 
that, if a A and an s P: \ were both determined, we could eventually use this 
information to recover the zero p. This might be used to write down a condition 
on the fixed points A that is equivalent to the Riemann hypothesis. The third 
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item we would need in this scheme is the function (p, A) cr p ,\ that takes a 
pair (p, A) to the arc-length of s p ,\ from A to p. We do not study cr p ,\ in this 
article. In fact no theorems are proved at all, but we do state conjectures based 
on computer experiments. The experiments were done with Mathematica and 
spot-checked with Sage. These experiments indicate that the backward orbits 
of the non-trivial zeros are controlled by a system of zeta cycles A. 

1.2 Definitions. 

Let ( denote the Riemann zeta function and let us write the iterates of a func- 
tion / as / o0 0) = z and f°( n+1 \z) = f(f on (z )) for n = 0,1,.... An n-cycle 
for / is an n-tuple (co, ...,c n _ i) such that /(c n _ i) = Co and f(ck) = c^+i when 
k ^ n — 1. The forward orbit of r under / is the sequence (r, /(r), / o2 (r), ...). 
The backward orbit of r under / is the set of complex numbers w such that 
f on (w ) = r for some integer n > 0. Let the symbol f°~(r) denote this back- 
ward orbit; if r does not belong to a cycle, f°~(r) carries the structure of a 
rooted tree in which the root is r and the children of w £ /° _ (r) are the solu- 
tions u of f(u) = w] obviously this is the case for f = ( and r = a nontrivial 
Riemann zero p. We will call any path in f°~(r) that begins at r a branch 
of f°~(r) (also: “a branch of the inverse of /”.) Such a branch, then, is a 
sequence (ao, ai, fl2, ••••) with ao = r and a n = /(a n+ i) for each non-negative 
integer n. Since the Riemann hypothesis has been verified within the range of 
our observations, we write without ambiguity pk for the k th nontrivial Riemann 
zero ordered by height above the real axis and p-k for its complex conjugate. 
If C on (z) = 0 and £ on_1 (z) is a nontrivial Riemann zero, we call 2 a nontrivial 
zero of £ on . For typographical reasons, we will occasionally write ( n for £ on ; 
within our article, there should no confusion with other common uses of this 
symbol. 

For 2 G C U { 00 }, A z := {u G C s.t. limn^oo ( on (u) = zj (the “basin of 
attraction” of 2 under zeta iteration.) Let 0 « —.295905 be the largest negative 
zeta fixed point. Then A $ is a fractal [11]; each nontrivial Riemann zero lies in 
an irregularly shaped bulb of A ^ (Figures 1.1, 3.2.) 

For a spiral s P: \ we define a real- valued function on complex numbers 6(z) by 
requiring that 9{z) = arg(z — A) (mod 27 r), and that 6{z) varies continuously 
and montonically as 2 moves around s P: \ in a fixed direction. In other words, 
0(z) behaves up to a multiplicative constant like a winding number. (For our 
immediate purpose, any such 6 will do.) We set r(z) = \z — A|. Suppose there 
are non-zero real numbers a, 6, c, and d such that all £ E s P: \ satisfy 

| log r(z) - (a + b0(z))\ = o(0(z)) 

with c exceedingly small; then we say that s P: \ is nearly logarithmic. If a 
sequence ( 21 , 22 ? •••) lies on a spiral s p , \ and for some pair of real numbers A > 
0, B > 1 and all k = 1, 2, ..., | arg(A — Zk) — arg(A — 2^+1 )| < Ae~ Bk , we will say 
that the Zk are distributed nearly uniformly around s P: 
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1.3 Summary of observations. 

We express some of our observations as explicit conjectures. 

Conjecture 1 For any positive integer L there is an infinite family Cl of Rie- 
mann zeta cycles A = (Ao, Ai, ..., Al-i) with the following properties. 

(1) For any non-trivial Riemann zero p and any A G Cl, each A G A is a 
repelling fixed point of £ oL lying in the boundary of A ^ under the usual topology 

on C. 

(2) Each A G A is the center of a nearly logarithmic spiral s p ,\, and the point 
set Uaga s pA contains a branch B p , \ of crip)- 

(3) The branch B p \ = (ao, ai, a 2 , ...) satisfies 

(1) a 0 = p, 

(ii) a n = C(a n+ i),n = 0, 1,2, ..., and 

(Hi) for each j = 0, 1,2, ...,L — 1, the sequences bj = (a^, a J+j L, Uj+ 2 L? •••) con- 
verge to A j. 

(4) The members of B p ^\ are distributed nearly uniformly on s P: \. 

(5) If z G C°~(p) then ( on (z ) G B p , a for some positive integer n. 

Because A j would be repelling, it would also be an attracting fixed point of a 
local branch of the functional inverse of and then the convergence would 
follow from standard results, for example, Theorem 2.6 of [4]. 

Conjecture 2 (1) The bj are interpolated by spirals \ with center A j. 

(2) The members of bj are distributed nearly uniformly around s^Xj. 

(3) The s p ^\ j are nearly logarithmic spirals. 

When L — 1, our experiments indicate that there is an infinite family of zeta 
1-cycles A = {A} of non-real repelling zeta fixed points A such that for any 
nontrivial Riemann zero p there is a branch B P: a = (ao, ai, a 2 , ...) of (°~(p) 
with ao = p, a& = C{ a k-\- 1) f° r all k > 0 and lima^ = A such that the a^ are 
interpolated by a nearly logarithmic spiral centered on A. If the forward zeta 
orbit of p, namely (p, 0, — 1/2, ...), is adjoined to B p 5 a to produce the doubly- 
infinite sequence 

E p \ — (• • • , a n , a n — i , . . .a \ , ao — p, 0, 2 ’ * •) 

then the entire ensemble E P: a appears to be interpolated by a double logarithmic 
spiral in the sense of [8], vol. 2, p.187, in which the center of the initial spiral is 
A and that of the terminal spiral is <fi. The point linking the two spirals, namely 
zero, is poorly approximated by both of them; the relative error is 200%. 

Conjecture 3 (1) There is an infinite set of repelling fixed points fi- 2 n near 
the trivial zeros —2 n < —18. 
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(2) The members of B p ^_ 2n y lie on a curve which is very nearly a straight line 

segment. If n is even , the endpoints are p and otherwise , the curve passes 

through ^- 2 n and terminates at p. 

(3) If2n = 0 (mod 4), then 


arg 


dQ^hn) 

dz 


27r; 


If 2n = 2 (mod A), then 


arg 


dCi^hn) 

dz 


(See Figure 5.7.) These observations are consistent, of course, with the propo- 
sition that B p ^^_ 2n y is interpolated by a spiral. Our computations seem to 
indicate that, while the imaginary part of ^- 2 n is always extremely small, it is 
not always zero. Conjecture 3 suggests that there might be a broader relation- 
ship of the same kind between the derivative of zeta at a given fixed point and 
the structure of the associated spirals centered at those fixed points elsewhere in 
the complex plane, but our experiments have not verified any such relationship. 
It is also suggestive, of course, that these relationships are exact when zeta is 
considered as a function of real numbers and ^- 2 n is replaced by —2 n. 


It seemed to us that there were two obvious possibilities to account for the 
observed connection between zeros and fixed points. It might exist, at bottom, 
because zeros of a function / are fixed points of s 1 — >• s + f(s) and fixed points 
of / are zeros of s H- f(s) — s. Or, on the other hand, the connection might 
not depend on special properties of the zeros at all. In other words, it might be 
the case that branches of (°~ (z),((z) ^ 0 converge to the same fixed points as 
branches of (°~(p),p a non-trivial zero. The second possibility seemed plausi- 
ble, as we will explain below, in view of an analogy between the spirals s and 
the streamlines of a fluid vortex. We tested various z on the critical line other 
than Riemann zeros and found spiral branches converging to the fixed points of 
zeta. So it ought be possible to explain the spirals with a theory that avoids any 
appeal to special properties of the Riemann zeros. We are not going to describe 
these experiments in any further detail in this article. 

The fluid analogy also led us to check for invariance of branches of (°~(z) under 
rotation about the fixed points at their centers. We found that the deviation 
from this sort of invariance is systematic and can itself be described by referring 
to (other) logarithmic spirals. 

We made a brief survey of functions other than zeta to gauge the extent of 
the spiral phenomenon, which we will not describe in any further detail than 
the following. Functions as simple as cosine appear to exhibit this behavior. We 
also observed it in, for example, the Ramanujan L-function. We hope to carry 
out another survey with a different software package. 
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1.4 Prior work. 

Many authors have examined the Riemann zeta function with computers. No- 
table citations from the perspective of this article are Arias- de- Reyna [3] , Broughan 
[1], Cloitre, [2], Kawahira [5], King [6, 7] and Woon [11]. 

2 METHODS 

2.1 Quadrant plots. 

We will be displaying colored plots (say, “quadrant plots” ) depicting, for a point 
w of C and a meromorphic function /, the quadrant of f(w). We use quadrant 
plots in three ways: (1) to determine small squares containing exactly one solu- 
tion of an equation of interest, so that this information can be used by standard 
equation-solving routines to find a solution to several hundred digits of preci- 
sion (which we find is necessary, for example, to locate zeta cycles) lying in a 
particular region; (2) to superimpose quadrant plots upon plots of the basin of 
attraction A^. These two kinds of plot typically interlock in a way that helps 
us to understand the meaning of many small irregular features of A and (3) 
to show how the quadrant plots spiral as we reduce the size of the plot window 
about a fixed point of zeta or one of its iterates. Observation of these spiral 
motions was our first indication that forward orbits near fixed points do lie on 
spirals. 

In quadrant plots, the boundaries of single-colored regions are / pre-images 
of the axes-curves corresponding to zero sets of $?(/(«§)) and ^(/(s)); the ap- 
parent intersections signal the presence of zeros or poles of /. By adjusting the 
color scheme to distinguish between regions where |/| is large or small, we can 
try to distinguish zeros from poles. Some apparent intersections are revealed 
to be illusions by a change of scale. Similar but colorless methods for plotting 
zeros were put to use in [3]. 

The visualized region is partitioned into small squares, each of which is rep- 
resented by a pixel. We choose a test point s in each square. The pixel rep- 
resenting the square is colored according to the rules in Table 1. In the table, 
the region D is a disk with center s = 0 and large radius r (chosen as may be 
convenient.) We denote the complement of D as —D. 
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Location of f(s) 

Color of pixel depicting region containing s 

real and imaginary axes 

black 

D n Quadrant I 

rich blue 

-Dfl Quadrant I 

pale blue 

D fl Quadrant II 

rich red 

— D fl Quadrant II 

pale red 

D fl Quadrant III 

rich yellow 

— D fl Quadrant III 

pale yellow 

D fl Quadrant IV 

rich green 

—D fl Quadrant IV 

pale green 


TABLE 1: COLORING SCHEME FOR QUADRANT PLOTS 


The junction of four rich colors represents a zero, the junction of four pale colors 
represents a pole, and the boundary of two appropriately-colored regions is an 
/ pre-image of an axis. An example is shown in Figure 2.1: (s- l) 2 (s — 

z)(s + 1 ) 5 / ( s + i ) 3 . (We have superimposed a pair of axes on this quadrant plot.) 



FIG. 2.1: QUADRANT PLOT OF 8 (s - l) 2 (s - i)(s + 1 ) 5 /(s + l) 3 
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2.2 High precision equation solving under geometric con- 
straints. 



FIG. 2.2: ZOOMING IN ON Ai 


We illustrate the application of quadrant plots to solve equations under ge- 
ometric constraints by showing how we found the three-cycle A described in 
section 5.4. In Figure 2.2, we have superimposed plots of A $ and quadrant 
plots of s i y ( o3 0) - son small squares near the first non-trivial Riemann zero. 
The resolution has been kept low to speed up the computations; high resolution 
is not particularly helpful in this situation. The left panel is a square with side 
length 20 and center pi. The central panel depicts a square with side length 2 
and center pi + 3.4 + .li; we have adjusted the center to keep in view a par- 
ticular four-color junction visible in the left panel. It represents a solution of 
( o3 (s) — 5 = 0, namely the three-cycle element Ai we are trying to compute. 
At this stage, if we used, for example, the Mathematica command FindRoot 
constrained to search within this square, it might land on any one of the several 
four-color junctions we see in the central panel. So we change the center of the 
plot again, this time to pi + 3.46 + .103i, and we make a square plot centered 
there with side length .02. This is shown in the right panel of Figure 2.2. Next 
we use a slow, “handmade” routine to approximate Ai by searching within this 
square. Then, using this approximation as the beginning value for a search with 
FindRoot, we obtain a solution with 500 digits of precision: 
3.9589623348847434673516458439896123461477039951866801455882506555054 
331235719797619129160432526832126428515417856326242408422124490775895 
37219976674458409141742662175701089081252727395073714398968532356378 
12255138302084634149524670809965144703541657360428502230820135428609 
38536894453944241116438492746243199878001238993540770158034816978947 
866042863811536518002674033394246742728451523022955079328623947833520 
567532298244004442294156837342370982002330874074322076777185746207730 
323482406094614280046 

+14.23622856322181332287122301085588169299871236208494399568695437825 

6069322896396761526007136189745757467102551375667154010366364994538731 

7916271113823253751110948972775762217941663830770714262041755664035323 

9671078789529204404394764315531582588051352309327292004343654135172820 





7780017861238006999109644383198471665302823015355865202971277187847669 

1974168218415293165267046606327405458655765280027732495125802150527924 

57282410834191507107658393848458313664113623935800293262678700791600125 

465010766853L 

Let us denote this approximation of Ai as a. A numerical check indicates that 
|C° 3 (a) — a | agrees with zero to 495 decimal places. 

Our main reason for requiring so much precision is that we will be repeatedly 
solving equations of the form ((u) = v for iq in each case replacing v with the 
previous u , to construct lists of (usually) 100 elements of a branch of the back- 
ward orbit of a nontrivial Riemann zero, looking for the ids near pre-selected 
fixed points. As the procedure repeats 100 times, there is an accumulation of 
numerical error, and in this situation very high precision is needed to maintain 
enough accuracy to “see” the spirals formed by these branches in our plots. 

3 A TOUR OF A 0 

We are interested in A $ because plots of this set make visible the underlying 
structure of the network of ( on pre-images of the critical line for all n at once: 
(1) the nontrivial zeros of the ( on he in bulbs of A ^ on filaments F decorating the 
border of A^, and (2) one ( on F pre-image of the critical line transects each such 
F. (Claims 1 and 2 are not, of course logically equivalent; we are summarizing 
computer observations that we will describe in more detail below.) Thus the 
structure of union of rooted trees visible in plots of A $ is apparently graph- 
isomorphic to a corresponding structure for the point set 

U C° _ (^) =W(say.) 

«(*)=§ 

This observation informs our discussion of the trees T in the next section. We 
pretend that we have stated a rigorous definition of the decoration notion and 
definite conditions for the membership of a given complex number in a given 
filament. In view of the relationship between U and A^, this should not cause 
problems: each filament F may be identified with one (of the many) ( ou f pre- 
images of the critical line, the definition of which could be made precise. But 
we should say explicitly that “A decorates F>” is a transitive relation and that 
the filaments are subsets of A<y 

In Figure 3.1, for example, the points at the junctions of four colors repre- 
sent zeros of £° 2 ; the zeros in the long filaments are nontrivial. The right panel 
of Figure 3.2 shows a quadrant plot of s i— ((s) — s superimposed on a plot of 
A^; the fixed points of zeta appear as the junction of four colors. The left panel 
depicts the nontrivial Riemann zeros using the same scheme (a quadrant plot 
of zeta.) 
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FIG. 3 . 1 : ZEROS OF £° 2 IN FILAMENTS DECORATING THE MAIN CARDIOID 



FIG. 3 . 2 : LEFT: ZETA ZEROS; RIGHT: ZETA FIXED POINTS 
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The filled Julia set of zeta (the points in C with bounded orbit under itera- 
tion by zeta) is C — A^. The basin A $ appears to be dense in C - A^. The 
sets C - Aqq (Figure 3.3) and A^, regarded as regions in the complex plane, are 
indistinguishable in our plots but they are not identical. For example, there is 



FIG. 3.3: C - Aoq 

an infinite number of real zeta fixed points ([11], Theorem 1) that belong to (C 
—Aqq) — A^. In addition, there appear to be infinite families of non-real zeta 
k - cycles for each integer k > 1 in (C —A^) — A^. 

Zero lies in A $ ([11], Theorem 1.) This set is a fractal decorated with numerous 
long filaments (Figure 3.1.) Zeroes of the ( on lie on the filaments. Because zero 
is an element of we know that the whole backward orbit £ o_ (0) H es i n A/>- 
Because the pre-images of nontrivial Riemann zeros under iterates of zeta lie on 
the filaments, the itinerary of a point in the backward orbit of a nontrivial zero 
p visits several filaments at the edge of A $ before coming to p. (Some but not 
all pre-images of the trivial zeros also lie on filaments.) 
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The set A $ seems to comprise 


(1) a heart-shaped, seven-lobed central body, which we will call the main car- 
dioid. 

(2) two major filaments of bulbs of various irregular shapes that emanate from 
the main cardioid, transected by the critical line and containing one nontrivial 
Riemann zero in each bulb (right panel of Figure 3.2.) 

(3) infinitely many blunt processes and long filaments decorating the main car- 
dioid and each of the irregular bulbs. The filaments comprise smaller copies of 
the bulbs, which, in turn, are decorated with similar filaments, ad infinitum. 
(Figure 3.1.) Thus, when we plot them, the set of filaments decorating A $ ex- 
hibit a visible tree structure. 

The visible features described in (1) - (3) were evident in Woon’s plots of C 
— A^ [11]. The filaments appear to be zeta-iterate pre-images (close copies) of 
the two major filaments. Pre-images of the real axis pass through the blunt pro- 
cesses and contain pre-images of the trivial zeros. For example, the right panel 
of Figure 3.1, superimposes a quadrant plot of ( ° 2 on the left panel, so that 
junctions of four differently-colored regions each represent a zero of ( ° 2 . There 
are three long filaments depicted in this image containing zeros, the immediate 
zeta images of which are nontrivial Riemann zeros; but between the lower two 
such filaments is a blunt process transected by a zeta pre-image of the real axis, 
and we can see another series of £ o2 -zeros lying along this curve. These are zeta 
pre-images of the negative even numbers. 

(4) at each trivial zero < —18, a microscopic, more or less distorted copy 
(zeta-iterate pre-image) of the entire assemblage described in (1) - (3). (By 
“microscopic” features we mean features so small that they can only be visual- 
ized by a change of scale from that of Figure 3.3.) In Figure 3.4, we show copies 
of the main cardioid near the trivial zeros —28, —26, —24, —22 superimposed on 
quadrant plots of £ ° 2 in the same squares. The size of these features decays 
exponentially with distance from zero. Their left-right orientation alternates. 
We speculate that the alternation can be derived from the alternating sign of 
the real derivative | a .__ 2n . 

Because these copies exist on the left half of the real axis, its zeta pre-images 
also contain complete copies of A^. The upper left panel of Figure 3.5 depicts 
the first bulb in the major filament in the upper half plane. It is a 10 by 10 
square centered at p\ . Along its border we see an apparently infinite set of fila- 
ments alternating with an apparently infinite set of blunt processes. Our tests 
demonstrate that the filaments are transected by £ on pre-images of the critical 
line for (n = 1, 2, 3, ...), and that the blunt processes are transected similarly by 
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fig. 3 . 4 : A ^ copies near s = - 28 , - 26 , - 24 , -22 superimposed on ( ° 2 

QUADRANT PLOTS 
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FIG. 3.5: A COPY OF A 0 NEAR A BLUNT PROCESS OF A 0 
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FIG. 3 . 6 : C U PRE-IMAGES OF THE CRITICAL LINE NEAR THE MAIN CARDIOID 
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( on pre-images of the real axis. The other three panels depict a small copy of 
A , ^ to the right of the largest blunt process on the right side of the bulb shown 
in the upper left panel. In the lower panels, a quadrant plot of ( 03 in the left 
panel and of £° 4 in the right panel have been superimposed upon this copy. 
Evidently, it is a ( 03 pre-image of A^. 

( 5 ) Our observations indicate that for each filament F decorating A ^ there 
is a positive integer kp (say, the degree of F) such that each bulb of F contains 
one nontrivial £ okp zero, and no nontrivial zeros of £ ok for any k ^ kp. Even 
under the Riemann hypothesis, it would not be necessary from first principles 
that degree k filaments are transected by £ o/c_1 pre-images of the critical line, 
even though that is the simplest possibility. But it seems to be the case. In Fig- 
ure 3 . 6 , the C o/c_1 pre-images of the critical line transecting degree k filaments 
decorating the main cardioid are shown for k = 1 , 2 , 3 and 4 . 

4 BRANCHES INTERPOLATED BY SPIRALS 

For any integer L > 0 there appears to be an infinite set of zeta cycles A = 
(Ao, ••., Al-i) that pick out a linearly ordered subset B P: a = (ao, ai, <12, •••) of 
C°~(p) such that (1) ao = p, (2) a n = £(a n +i),n = 0,1,2,.... and ( 3 ) for each 
j = 0, 1,2, ...,P — 1, the subsequences bj = (aj, a J+J L, %+2L, •••) appear to con- 
verge to A j. The sequence bj is a branch of CL _ ( a i)* f n mos t °f the cases that 
we have examined, each bj appears to be interpolated by a spiral s pNi \ j with 
center Aj = lim bj. The A j are repelling fixed points of £l- 

Now we offer (speaking loosely) a geometric description of some of the bj in 
terms of the basin of attraction A^. (It applies to most, but not all, instances 
we have examined to date.) A variety of filaments decorate A^, but here we 
restrict attention to those that decorate the main cardioid. We assign the set 
of filaments a structure of union of rooted trees T as follows. A filament FgT 
is the parent of a filament G G T if and only if G decorates F and there is no 
intermediate filament H E T such that G decorates H and H decorates F. The 
filaments containing nontrivial Riemann zeros have no ancestors, but they are 
not unique in this respect. 

Now fix integers L > 1 , 7 V 7^ 0 . There is an infinite set of zeta P-cycles 
A = (Ao, A2, ••., Al-i) such that for each integer A = 0, 1,2, ...,P — 1, there 
is a tree Xa of filaments decorating the main cardioid of A^, and a path 
Pa = (Fo, Pi ,-•...) in Pa with kp m = A + mL and such that (if m > 0 ) F m 
decorates the | 7 V| th bulb of its parent F m _ 1. As in the first column of Figure 
4 . 1 , the filaments in Pa spiral around Aa- In our graphic visualizations, the 
apparent size of F m decays exponentially with m. Something like this would 
seem to be a necessary condition of the relation A j = lim bj we mentioned above. 

Each bulb of F m contains a nontrivial zero w of £ oA + mL and ( oA ^ rnL ~ 1 (w) 
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FIG. 4 . 1 : FILAMENTS SPIRALING AROUND 1p pi 
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is a nontrivial Riemann zero. Which one? Let re m;j /v be the nontrivial £ oA + mL 
zero belonging to the \N\ th bulb of the filament F m in Fa- For m > 0, 
so the sequence (wo,n, ...) is a branch of CiT( w Q,n)- 
In our observations, the zeta image of bulb \N\ of a filament F in Ta with kp > 1 
is bulb | TV | of its parent filament in Ta- Therefore C° A ^ rnL ~ 1 (wm, N ) = p±N- 

The observation that = w m - i,iv is our evidence for an idea men- 

tioned in the introduction, namely, that there should graph isomorphisms be- 
tween subgraphs of the rooted tree graphs associated to the on one side 
and subgraphs of the trees T decorating the main cardioid of A ^ on the other. 
We should mention that the situation for copies of the main cardioid such as 
the ones illustrated in Figures 3.4 and 3.5 is different; except to say that the 
zeta images of copies are also copies, we will not discuss it further in the present 
article. 

5 SPIRALS INTERPOLATING A BRANCH OF C~(p) 

5.1 Single spirals interpolating a branch. 

When L = l, A = {Ai} where Ai is a repelling zeta fixed point; there appear 
to be at least three categories of such points: ^- 2 n (say) lying near the trivial 
zeros —2 n = — 20, — 22, ...; zeta fixed points ip p * near each nontrivial Riemann 
zero p* , and eight fixed points lying at the boundary of the main cardioid (right 
panel, Figure 3.2.) How near? In the case of the p*, one can form an impression 
by keeping in mind that this figure depicts a 120 by 120 square (section 8.) The 
distances | — 2 n — ^- 2 n\ are a great dealer smaller; we omit the details. 

There are exactly two filaments F with degree kp = 1 decorating the main 
cardioid; one of them contains pN in its \N\ th bulb = /3/v, say. Our observations 
are consistent with the following proposition. The point 'ip p N lies at the border 
of the | TV | th bulb of a filament F* with kp* = 2 decorating /3/y This ramifies: 
if vp PN lies at the border of the \N\ th bulb of a filament F* then there is a child 
filament F' of F* such that ip PN lies at the border of the |7V| th bulb of F' and 
kp' = kp* + 1. 


This is illustrated by the left column of Figure 5.1 for N = 1, 2, 3, 4. It depicts 
quadrant plots of s i— £(s) — s, so that r ip PN shows up as four-color junctions 
on the depicted filament (say, F 2 .) The right column shows quadrant plots of 
C o3 ,C o4 ,C o5 ,C° 6 i n rows 1, 2, 3 and 4, respectively, all superposed on plots of 
Afj,. The squares have side length .2, .02, .002, and .0002 in rows 1, 2, 3 and 
4, respectively. The center of the squares in row N is 7 /^, so the panels are 
depicting the region around this point at smaller and smaller scales. 

Figure 4.1 displays the original indications we had that some branches of the 
inverse of zeta lie on spirals. It zooms in on the illustration of ^p pi in the left 
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FIG. 5.1: Wn> 1 < n < 4 


19 





panel of the top row of Figure 5.1. The center of that panel is the fixed point, 
lying on the border of a filament F% (say) decorating the lower border of the 
largest full bulb. In the right panel of the top row of Figure 4.1, a quadrant plot 
of C° 3 has been superimposed on a plot of the corresponding region of A we 
see that F 3 contains zeros of (° 3 . (Simple tests show that they are nontrivial 
zeros in the sense of the introduction.) In the lower rows, the zoom is repeated 
and VVi i s seen t° he near a still-smaller filament decorating the bulb near the 
center of the figure just above it. The right column depicts quadrant plots of (° 4 , 
c 5 and ( 0 6 for the squares opposite them in the left column. So in Figure 4.1 
we are seeing zeros of these functions (again nontrivial.) They also appear (in 
virtue of the shapes of the underlying A^-bulbs) to be C o3 :C° 4 :C° 5 pre-images 
of nontrivial Riemann zeros. The rapid reduction of scale from one row to the 
next attests to a similar reduction of the distances of these pre-images from ^ pi 
(which, as we have remarked, is not surprising.) The possibility that they may 
be traveling on spirals emerges from a look at the angles that the filaments F% 
and (say) F^F$,Fq make with the horizontal. These observations led us to do 
the numerical tests described in the last section. 

We made a survey of the spirals s p ^ p * for various choices of p and p*. We 
made a table of ijj Pn , 1 < n < 100 with 500 digits of precision and we used a 
table of nontrivial Riemann zeros with 300 digits of precision made by Andrew 
Odlyzko [9]. We made tables of the zj~ in various to high precision, 

proceeding inductively. We set zi = p and, given a value of Zk, after using 
Mathematics s FindRoot command to solve ((s) = Zk in the vicinity of Vy>*, we 
set Zk+i equal to the solution. We began with 300 Zk for each B and used tests 
of reliability of each Zk to truncate the list; typically, we ended up with a least 
100 consecutive z\~. 

We use polar coordinates (r(z),0(z)) to denote a typical point 2 on <s p ,p p * 
such that (1) r(z) = \z — Vvl, (2) 0 ( z ) is chosen so that 6 (z) = arg(z — ^ p *) 
(mod 27r), and (3) 6 (z) varies continuously and monotonically as z moves around 
the spiral in a fixed direction. In other words, 6 (z) behaves up to a mutiplica- 
tive constant like a winding number. Then r(z) appears to decay exponentially 
with 0(z). (Of course we are only able to check this for 2 G F?p,p p * , that is, for 
the Zk we propose are interpolated by s p ,p p * , because no other test for mem- 
bership in Sprf * is available to us.) Therefore it was not practical to plot the 
spirals s p ^ + without re-scaling, so we plotted the points (logr(z), 0 (z)) instead. 
This procedure everts the apparent spirals: if k,j are such that r(zk) < 1 and 
r(zj) < 1, then r(zj) < r(zk) implies that the plotted point (log r(zj), 6 (zj)) 
is further from the center of the re-scaled interpolating spiral than the point 
(log r(zk),0(zk)): the reverse of the situation before re-scaling. (The direction 
of winding of the spiral is also reversed because the logarithms take negative 
values.) The points near the center of the re-scaled spiral depict Zk for smaller 
values of k for which Zk is closer to p and further from p* . They are crowded so 
closely, in spite of our re-scaling, that the interpolating curve near p is obscured. 
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Figure 5.2 depicts branches of backward orbits of p n (l < n < 5) spiraling 



FIG. 5.2: BACKWARD ORBIT BRANCHES FOR p n (1 < n < 5) CENTERED ON 
TWO CARDIOID ZETA FIXED POINTS 

around two fixed point « —14.613 + 3.10H (left column) and « —5.28 + 8.803i 
(right column) on the border of the main cardioid; we omit the 500 digit deci- 
mal expansions, which are easy to compute using Mathematical s FindRoot com- 
mand. Figure 5.3 depicts branches of backward orbits of p n spiraling around 
VV n (1 < n < 10). (We omit their precise expansions for the same reason.) 


5.2 An example: the case p — p* = p\. 

The first panel of Figure 5.4 plots the point set B pu ^ (re-scaled as described, 
and shifted to place the apparent spiral’s center at the origin.) It offers the 
appearance that the z k (in red) form arms something like those of a spiral 
galaxy; this seems to be a result of nearly regular growth of 0(zk) with k. But 
Zk for consecutive k do not lie in adjacent positions on these arms; in the second 
panel, the Zk are connected by chords in the same order as they appear in the 
sequence vertices v,w representing z v ,w = (( z v ) are connected by a 
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FIG. 5 . 3 : BACKWARD ORBIT BRANCHES FOR p n NEAR 1 < n < 10 



FIG. 5 . 4 : LOGARITHMIC SCALING OF s pi , pi 
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chord. The lower left panel of Figure 5.4 is a plot of log | z k — VVi | vs. it is 
clear that the distances of the z k (colored red) from the fixed point ip Pl at the 
center of the spiral are decaying exponentially. The lower right panel of Figure 
5.4 depicts a spiral curve (colored blue) that approximately interpolates the z k \ 
we found this curve using the NonLinearModelFit command in Mathematica. 
The equation of the curve is 

log r(z) = a + bO(z) + cexp(dO(z)), 


with 

a « 0.05575203301551956560399459579161529353, 
b « -2.39481894384498740085074310912697832305, 
c « -2.8680355917721941635331399485184884 x 10“ 120 , 

and 

d « 0.97375124237020440301256901292961731822. 

The parameter c is so small that this is quite close to being the equation of a 
logarithmic spiral. In the section on error terms below we will compare directly 
the loci of the z k with logarithmic spirals. 


In Figure 5.5, we study the variation in 6(z &). The left panel plots 4 = 
6(zk+i) — 0(z k ) against k and shows that the 0(z k ) are very nearly periodic 
in k. The right panel, which plots log |4+i — 4| against &, shows that the de- 
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FIG. 5.5: LEFT: S k VS. k RIGHT: log |4+1 -41 VS. k. 


parture from periodicity in 0(z k ) actually appears to decay exponentially with 
k. However for other choices of p and p* this no longer holds, and so it is an open 
question whether or not it would hold even in this example for very large &, that 
is, very close to the center of the spiral. We remark that the z k could, of course, 
be distributed along a nearly-logarithmic spiral while also being distributed in 
a completely irregular or at least non-periodic way in the theta aspect, so the 
two questions are at least superficially independent. 


Now suppose (ri, 0i) and (7*2, 4 ) He on a true logarithmic spiral logr = a + bO. 
The constants a, b are determined by any two points of the spiral, hence, if two 
pairs of points determine different values of a and 5, then the curve that the 
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three (or four) points comprising the pairs he on is not a logarithmic spiral. We 
used this idea to test B pi ^ pi = {21,22,23? •••} for the property of being inter- 
polated by a logarithmic spiral. We performed the test by solving for a and b 
using the pairs (21, z k ), k = 2, 3, .... In Figure 5.6 we have plotted the resulting 



FIG. 5.6: PARAMETERS a, b FOR LOGARITHMIC CURVE log r = a + bO INDUCED 
BY SUCCESSIVE POINTS OF B pupi . LEFT: a VS. b ; CENTER: b VS. fc; RIGHT: 

a vs. k 


values of a against 5, b against k , and a against k. Evidently a is roughly linear 
in 6, and both a and b appear to converge as k grows without bound. Thus 
the interpolating curve is not a logarithmic spiral, for then a and b would be 
constants. But the convergence of a and b suggests that the interpolating spiral 
s P i,^ P1 resembles a logarithmic spiral more and more closely as it winds inward 
towards x/j pi . 


5.3 Backward orbits near the trivial zeros. 

There appear to be zeta fixed points ^-2 n near each trivial zero —2 n < —18. 
Whether they lie slightly to the right or to the left of —2 n along the real axis 
appears to depend upon the parity of n. This reflects the alternating left-right 
orientations of copies (zeta pre-images) of the basin of attraction A $ we see in 
Figure 3.4. As we have already remarked the imaginary parts of the ^-2 n are 
small, but they probably are not all zero. 

A branch B p ^ of the backward orbit of each nontrivial Riemann zero p lies 
on a curve appearing to pass through or terminate at ^-2 n' if — 2n = 0 (mod 
4), then the curve appears to terminate at ^-27 These curves closely resemble 
straight line segments. Error terms are discussed in section 6. If — 2n = 2 (mod 
4), then (supposing, for the moment, that the curve really is a line segment) 
r ip— 2 n lies near its midpoint. In Figure 5.7, several of the backward orbits are 
depicted, re-scaled logarithmically as above. 

This observation is consistent with the hypothesis that B p ^_ 2n is interpolated by 
a spiral such that the a k G B p ^_ 2n satisfy | axg(a k -^-2n) -arg(a/c + i -^-2n)| « 
27 r for all k if —2 n = 0 (mod 4), or | arg (a k — ^-2 n) ~ arg(R/c+i — ^- 2 n)\ ~ 7r if 
— 2n = 2 (mod 4). We discuss this further in sections 5.5 and 6.2. 
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FIG. 5.7: COLUMN 1: BACKWARD ORBITS OF p 1 NEAR ^-2n, 10 < n < 13 
COLUMN 2: BACKWARD ORBITS OF p n _g NEAR ^_2 n, 10 < 7i < 13 
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FIG. 5 . 8 : QUADRANT PLOTS OF S ^ ( on (s) - S, 2 < n < 5 , NEAR THE BULB 

AROUND pi 
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5.4 Several spirals that together interpolate a branch. 

Figure 5.8 depicts members of zeta n-cycles (zeros of s i-A C° n ( 5 ) — s ) near the 
bulb of containing p% for 2 < n < 5. As n increases the pattern of distri- 
bution of these zeros becomes more and more obscure. The situation near pi 
appears to be typical of that near all nontrivial zeros of zeta iterates. 

Figure 5.9 illustrates the branch B Pi ? a of the backward orbit of (°~ ( pi ) induced 



FIG. 5.9: THE BRANCH OF A pi BACKWARD ORBIT INDUCED BY A ZETA 

3- CYCLE 

by a 3-cycle A = (Ai, A 2 , A 3 ) with Ai ~ 3.95896 + 24.2362i . The red vertices of 
a given chord of the graph represent points in the branch. Geometrically, Figure 
5.9 is doubly abstract: (1) The spirals have been positioned so that their centers 
are placed at 1000 + 0 i, 1000 (cos(^ L ) + i sin(^ L )) and 1000 (cos(^ L ) + i sin(^ E )) 
for the sake of legibility; and ( 2 ) the spirals are everted: they have been re-scaled 
logarithmically, so that points closer to the center appear, in the figure, to be 
farther from the center. The row 2, column 2 panel shows vertices representing 
elements of F> Pi ,a and edges between vertex pairs (v, ((v)), while the other pan- 
els depict the spirals separately; these are portraits of bj(j = 1,2,3). In these 
three figures, each vertex pair (v, ( 03 (v )) is connected by an edge. 

5.5 Angular distribution of branches along the spirals. 

A structural invariant, which seems to determine the number of arms visible 
in our plots of branches of /° _ (z), is the function Sf :Z ^ : k arg(a/ c — - 0 ) — 
arg(a/ c+ i — - 0 ), where the a/ds are members of a particular branch of f°~(z ) 
converging to ijj. 
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The values of 5k in the case discussed in section 5.2 correspond in the present 
notation to those of 5^ Pl? (k). They appear to change very slowly with &, and 
this behavior seems to be what gives rise to the appearance of discrete arms in 
plots of branches of C°~{pn )• 

Among the zeta fixed points very near trivial zeros, the only values of 
^C,pn,b_ 2 n— is (^) that we see (Figure 5.7) are ~ 7 r and ~ 27 t, distributed, as we 
have noted above, according to the mod 4 residue classes of the zeros. 

Because 5^ Pn ^ pn (k) apparently converges rapidly as k oc (as in Figure 5.5 

where n = 1), we take the value of 5^ Pn ^ pn (100) as a proxy for 

lim/^oo 5^ Pn ^ pn (k). Then our calculations are consistent with the proposition 
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FIG. 5.10: 
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1 < n < 600 


that lim/^oo S^ iPnj7 p ~ tt/2 (Figure 5.10.) Very small differences in this limit 
as n varies appear to determine very different shapes for the discrete arms visi- 
ble in our plots. 

We have observed in all of our experiments that the visible structure of a branch 
of f°~(z) depends upon the fixed point at its center and not on z, so 
should depend only upon / and i/j. Contrary to the impression suggested by our 
notation, it should be independent of z, but we cannot exclude the possibility 
that there are counterexamples to this idea. 
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5.6 Logarithmic models of spirals interpolating branches 
of the backward orbit of zeta. 

The branches B p ^ = (ao = p, ai,a 2 , ...) of (°~(p) f° r nontrivial Riemann zeros 
p converging to zeta fixed points V’ are interpolated by curves that resemble 
logarithmic spirals. We carried out experiments in which we looked for approx- 
imations of these interpolating curves by such spirals. We chose branches of 
the argument function, varying with k and evaluated at a — i/j, such that an 
angle Ok was assigned to ^ which was the least such angle > c + max J< / c Oj 
for c = 0 or 1. The angle Oo was the value of arg(ao — VO from the branch of 
argument chosen automatically by Mathematica. The choice of c was dictated 
by requiring that Ok act like a winding number about evaluated at the points 
ak . For V> near a trivial zero, cm 1 was chosen; for V> near a nontrivial zero, 
c = 0. 

For rk = \dk — V’l plots of the sets of pairs (0fc,logrfc),fe > 0, appear to lie 
on curves resembling straight lines. For spirals centered at zeta fixed points 
Vy> n near the p n , we approximated these lines using Mathematical s FindFit 
command. Figure 5.11 is a pair of plots of m n and b n against n for models 



fig. 5.11: slopes m n (left panel) and intercepts b n (right panel) in 

LOG-LINEAR MODELS T = exp(m n 6> + b n ) OF SPIRALS INTERPOLATING 
BRANCHES OF p n CENTERED AT 'lfj pn WITH r = \z — ^ Pn | AND 0 = arg(^ - ^ Pn ) 

l z ~ = oxp (m n 0 + b n ) fitted to branches of C° _ (Pn) 5 1 < n < 600. In par- 

ticular, we write fk = exp (m n 0k + b n ) for our estimate of These seem to be 
first-order approximations to genuine interpolating curves; the error-term will 
be discussed in the next section. 

Investigating spirals centered at the zeta fixed points ^ 2 n that lie near the 
trivial zeros was carried out in a different way. We were able to collect a sub- 
stantial amount of data (meaning for the first 200 members of the branches) 
for the VVn ^ n < 600 using 500-digit precision. On the other hand, even with 
1000-digit precision, we were able to collect data only on the first twenty ele- 
ments of branches centered at the V>-2 n fc> r n, < 30 before the use of the FindFit 
command to get a linear model for the pairs (0k, logr^) produced error messages 
from Mathematica. 
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Fortunately, the branches spiraling about the ^- 2 n appear to be better be- 
haved than the ones spiraling about the ip Pn • For each pair (n, n*) the branch 
B Pri ^_ 2ri * is apparently interpolated both by a curve very nearly a logarithmic 
spiral and by another curve which is very nearly a straight line passing through 
the points p n and i/^n* • As we will see in the next section, the fit of the 
branches to the straight line passing through these two points is so good that 
we use it as our second model, together with the assumption that the 0 k = cork + 
(a constant depending only on n and n*), with a — 1 or 2 depending as we have 
explained only on the parity of n*. It was feasible to find linear models for 
the maps k log r k - Combining the assumptions about the Ok with the linear 
models we construct for the log r k gives logarithmic models for the interpolating 
spirals. We test these models in the next section. 

We chose to examine the behavior of branches B Pn ^_ 2n _ 18 of C°~(pn) because, 
among zeta fixed points close to the trivial zeros —2 n, the greatest one (the one 
that lies rightmost along the real axis) is very close to —20. Figure 5.12 is a plot 
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FIG. 5.12: slopes m n (left panel) and intercepts b n (right panel) 

IN MODELS fk = exp (m n k + b n ) OF SPIRALS INTERPOLATING BRANCHES OF p n 
CENTERED AT ^-2n-18 

corresponding to Figure 5.11 for the zeta fixed points '0_ 2n _i8 . Here we have 
B Pn , , 0_2n— is = (ao = Pn,Q‘i,a 2 , •••)• Writing r k = |^- 2 n-i 8 - a k \, we took m n 
and b n to be, respectively, the means of the slopes and intercepts of the chords 
connecting consecutive pairs p k = (k,\og\a k — V ; - 2 n-i 8 |)* Thus our model for 
r k = \cL k — , 0_2n— 18 1 is fk = exp (m n k + b n ). We discuss it further this in the 
following section. 
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6 ERROR TERMS 

We will take the phrase “error term” to encompass complex-valued deviations 
from a given estimate as well as their absolute values. Like the original estimates, 
the curves followed by complex-valued deviations appear to have the form of 
logarithmic spirals. This raises the prospect of an infinite regress, which might 
perhaps lead to an exact expression for the best interpolating curves, but we 
have postponed any investigation of this idea. 
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6.1 Deviation of backward orbits from logarithmic spirals. 

6.1.1 Branches converging to fixed points near non-trivial zeros. 

We plotted the relative complex- valued deviations d(n, k) := 

(ak -t/yj -f k e iek 

a k — V’ p„ 

of the members of B p ^ from the logarithmic spirals approximating the best 
curves interpolating B p The runs depicted in Figure 6.1 portray both kinds 



FIG. 6.1: B Pri ^ n ,n = 1,28,48: COLUMN 1: ORIGINAL BRANCHES; COLUMN 2: 
DEVIATIONS OF B pn ^ n FROM LOGARITHMIC SPIRALS; 

COLUMN 3: log d(n,k) VS. k 

of plots for p = p n , -0 = ^ Pn for n = 1, 28, and 48. It seems noteworthy that the 
two kinds of plots resemble each other so closely, but inspection demonstrates 
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that they are not identical. 


The values of log d(n,k) for n = 1,28, and 48 are also plotted in Figure 6.1 
(column 3.) The magnitude of the d(n, k) appears to decay exponentially for 
k < roughly 130; for larger &, however, the magnitude of the deviations appears 
to grow exponentially without exceeding e~ 6 for k < 200. We think that the 
shape of this curve, which is typical, is an artificial effect of the FindFit com- 
mand on a file of 200 pieces of data: the fit is best near the center of the data file. 

For a fixed point = ijj Pn near a nontrivial zero p = p n , we used the initial 200 
elements of each branch B p as a proxy for B p to study the relative deviations 
of the curve interpolating it from a logarithmic spiral. Taking /? = 200 for the 
moment, let us set 

max nJ 3 := maxi < k <pd(n,k), 

max n,p := A/n/logn x x n ^, 

let mean n ^ denote the mean of d(n, 1), d(n, 2), ..., d(n, /?), and let 
mean* n p := sj n/ log n x \i n #. 

The panel in row 1 column 1 of Figure 6.2 is a plot of 600 values of mean* n 



FIG. 6.2: mean* o (left) and max ^ n (right); (3 = 200, 1 <n< 600 

The panel in row 1 column 2 is a plot of max* The plots in row 2 are smooth- 
ings of the plots in row 1: for each n, they depict means of mean*p and mar*^ 
over the range 1 < j < n. 
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These plots are consistent with the proposition that, for /3 = 200, fi n ^ and 
both = O with both implied constants < 1. More optimistically, 

perhaps, the plots are consistent with the hypothesis that mean n ^ and max n ^ 
both = o • Referring to column 3 of Figure 6.1, it seems possible that 

for fixed n the limits as (3 — > oo of max n $ and p n ,p both exist and that both 
= O or o as we ^* One reason to believe this is that in other 

ways the tails of the sequences B p ^ behave reasonably, as in Figures 5.4 and 5.5. 


We tested the same idea after replacing with powers for \ < 

e < 1. It seems possible that the supremum of e for which these statements 
might be true lies in the half-open interval [.8, .9). It also seems possible that 
this supremum is a decreasing function of n. We omit the relevant plots. 


6.1.2 Branches converging to fixed points near the trivial zeros. 

Let the branch B Pn ^_ 2n _ 18 of C°~(pn ) = (&o = pn,a i,a 2 ,....). Before we test 
a logarithmic model for the decay of = |a& — i/j- 2 n -is|, we want to assess 
how well B Pn ^_ 2n _ 18 fits the straight line passing through p n and ^_ 2n -i8- We 
measured the vertical deviation of the G B Pn ^_ 2n _ 18 from the straight line 
passing through both p n and ^- 277 - 18 as a fraction of the heights of the a^. We 
make the following definitions: 

M n and B n are the slope and intercept respectively of the straight line passing 
through p n and ^- 271-18 


and 


n, k) := 


- (M n 5ft(qfc) + B n ) 
3 ( a k ) 


mean t ™£' Lal = the mean of the d trivial (n , fc), 1 < k < /?, 
and 

max^ ial = maxi<K^ d trivial (n , k). 

The left panel of Figure 6.3 is a plot of log mean t r [' L p ' ial , 1 < n < 32 and (3 = 20. 
The right panel is a corresponding of log max l ^ ial . Evidently, the a & lie near 
the specified lines, and agreement with the lines improves rapidly as n increases. 


33 



-40 


-60 

-HO 


-100 



-60 


-100 



• -12 0l .... . ♦ 

5 10 15 20 25 ’ 30* J 5 10 13 20 25 30 * 

FIG. 6.3: log mean l ^ ial (left) and log max l ^ ial (right); 

’ p = 20, 1 < n < 32 


Next we define for a k e B Pn ^_ 2n _ 18 

^n,k — |R/c 2n— 18 1 5 

for k > 1, m n? / e = the slope of the chord connecting the ordered pairs (&, log r n , k ) 
and (&, log r n ,k-i) in R 2 , 

m n ^p = mean of the m n? / c , 1 < k < /3, 


and we define a ^/-intercept function b n ^ analogously. The error functions are 
defined as follows: 


d model {n,k,/3) 


logF n ,/c - {rrin^k + b n ,p) 
log r n ,k 


mean™°p del = the mean of the d modeZ (n, fc, /?), 1 < k < /?, 
and 

max™ del = maxKK^ d model (n, fc, /3). 

The left panel of Figure 6.4 is a plot of log mean™ odel , 1 < n < 32 and /3 = 20. 
The right panel is a corresponding of log max ™ odel . Once more the fit is good 
and improves rapidly as n increases. 

7 An analogy from fluid mechanics. 

It did not seem plausible to us that something special about Riemann zeros p 
should force the branches £° _ ( 2 :), z = p in particular, to be attracted to repelling 
fixed points ip along logarithmic spirals. As we remarked in the introduction, 
dynamical systems theory leads one to expect that the pj should attract all of the 
nearby branches (°~(z) whether or not ((z) = 0, and (one speculated) probably 
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FIG. 6.4: log mean™°p del (left) and log max™°p del (right); 

’/3 = 20,1 < n < 32 



along roughly similar curves. Logarithmic spirals appear in fluid mechanics (see, 
e.g ., [8], v.2, pp. 186-188 or [10], p. 358.) By analogy with the streamlines of 
a vortex in a fluid, we speculated that the existence of spiral curves connecting 
zeros of zeta to repelling zeta fixed points might be a consequence of a scenario 
in which there is an infinite family of such spirals related by rotations around 
the fixed point. By this we mean a family of spirals s parameterized by real 
numbers x, varying continuously with x in a sense made explicit by condition 
(1) below, such that if s x and s x +q are two such spirals with common center a 
zeta fixed point ip, then 


(1) s x+e -ip = e i0 (s x -ip) 


(condition (1) being an equation of homotheties) and 

(2) z E s x => (i) ((z) E s x and (ii) there exists a branch B z C s x of (°~(z) such 
that lim B z = ip. 

In this scenario, the spiral curves would be congruent in the sense of Euclidean 
geometry and exactly one spiral would intersect the critical line at each p n with- 
out appealing to special properties of the zeros. We have verified the existence 
of spiral branches B z of (°~(z) for various z on the critical line other than 
Riemann zeros without meeting a counterexample. Like the spirals we have al- 
ready described, they are approximately logarithmic; we omit the relevant plots. 

Condition (1) imposes rotational invariance on the s x . This suggests the possi- 
bility that the branches of (°~ interpolated by them enjoy the same property. 
Suppose u and ((u) lie on s x with center ip and let Rq^(z) := e' l0 (z — ip) + ip be 
the function that takes 2 to its image under rotation by an angle 0 around ip • 
Under perfect rotational invariance, not only of the spirals s x but of the branches 
of (°~(z) for particular z that they interpolate, the numbers Rq^(((u)) — 
((Ro^(u)) must vanish. Therefore we studied the Re^(((u)) — ((Rq^(u)). 
We restricted ourselves to u E C°~( z ) f° r various z, not only because these are 
the main objects of interest, but because our only reliable information about 
the spirals comes from their interpolation of the branches B ^ = (ao, ai, e&2, •••) 
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of (° (z), and so our only useful candidates for points in s are the members of 
such branches. 

Logarithmically scaled plots (which we omit) of the discrepancies (say) Ro^(C(a n )) — 
C{Re^(a n )) indicate that these numbers decay in modulus exponentially and 
rotate around the origin in a nearly linear fashion with n. In other words, they 
themselves describe curves that are approximated by logarithmic spirals. 

8 APPENDIX: THE FIGURES 

8.1 Figure 1.1 

Figure 1.1 a 120 by 120 square with center 1 + Oi. 

8.2 Figure 2.1 

Figure 2.1 depicts a 6 by 6 square with center zero. 

8.3 Figure 3.1 

Figure 3.1 shows an 8 by 8 square with center — 5 + 9.5i. 

8.4 Figure 3.2 

Figure 3.2 shows a 120 by 120 square with center zero. 

8.5 Figure 3.3 

Figure 3.3 depicts a 60 by 60 square with center zero. 

8.6 Figure 3.4 

The upper left panel of Figure 3.4 shows a 2.4 x 10 -5 by 2.4 x 10 -4 square 
centered at —28. The upper right panel shows a 2.4 x 10 -4 by 2.4 x 10 -4 square 
centered at —26. The lower left panel shows a .004 by .004 square centered at 
—24. The lower right depicts a .07 by .07 square centered at —22. 

8.7 Figure 3.5 

The panel in row 1 column 1 of Figure 3.5 depicts a 10 by 10 square centered 
at p\. The other panels show a square with side length .006 and center pi + 
4.1215 - .4015i « 4.6215 + 13.7332*. 

8.8 Figure 3.6 

Each panel of Figure 3.6 is a 30 by 30 square with center = — 5. 
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8.9 Figure 4.1 

The squares depicted in Figure 4.1 have side length .2, .02, .002, and .0002 in rows 
1, 2, 3 and 4, respectively. The center of each square is 'ip pi ~ —2.3859 + 16. 271b 

8.10 Figure 5.1 

All the panels of Figure 5.1 depict A ^ in 2 by 2 squares. In rows 1-4, the 
centers are x/j pi - ^ P4 , respectively, where Vb 2 ~ —2.0369 + 21.993H, !p P3 ~ 
-1.6935 + 26.5283i, P and xp P4 ~ -1-7496 + 30.8158b 

8.11 Figure 5.2 

In Figure 5.2, column 1 depicts branches of C° (Pn), 1 < n < 5, centered at a 
zeta fixed point « —14.613 + 3.108i; column 2 depicts branches of C° _ (/+)> 1 < 
n < 5, centered at a zeta fixed point ~ — 5.279 + 8.803b 

8.12 Figure 5.3 

In Figure 5.3 (referring to the caption), the value of n in row a, column b is 
n = 2n + 6 — 2. 

8.13 Figure 5.8 

Figure 5.8 depicts four views of a 12 by 12 square with center pi. 

8.14 Figure 5.9 

Figure 5.9 depicts the branch B p 5 a of the backward orbit of C °~{pi ) induced 
by a 3-cycle A = (Ai, A2, A3) with Ai « 3.95896 + 24.2362b 
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